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Abstract. We construct a new symplectic, bi-hamiltonian realization of the 
KM-system by reducing the corresponding one for the Toda lattice. The bi- 
hamiltonian pair is constructed using a reduction theorem of Fernandes and 
Vanhaecke. In this paper we also review the important work of Moser on the 
Toda and KM-systems. 



1. Introduction 

Some of the most important work of Moser concerns the Toda and Volterra lat- 
tices. These two systems are closely related and Moser gave an explicit construction 
demonstrating the relationship between the two systems. The Toda lattice was dis- 
covered by Morikasu Toda in 1967. Computer experiments by Ford et al. [22] 
suggested that the Toda lattice is integrable. In two papers, both in the same issue 
of Physical Review B ([19J, |23J) Henon and Flaschka in 1974 demonstrated the 
integrability of the system. Henon provided the required independent constants of 
motion using combinatorial arguments. Flaschka proved integrability of the lattice 
via a change of variables and the construction of a Lax pair. A few months later, 
Manakov [27] established independently similar results. The Toda lattice is a dis- 
crete approximation of the KdV equation. In [2U] Flaschka used this analogy to 
solve the system using a discretization of the inverse scattering method of Gardner, 
Greene, Kruskal and Miura. The method which was developed for partial differen- 
tial equations was applied for the first time for the infinite Toda lattice with success. 
Moser's paper [28] deals with the finite non-periodic Toda lattice. Moser solves the 
system by constructing action-angle coordinates. His action variables are the eigen- 
values of the Lax matrix of Flaschka and the angles are the last components of the 
normalized eigenvectors. We will describe the solution of the Toda lattice by Moser 
in Section 2. The Volterra system (also known as the KM system) is closely related 
to the Toda lattice, and the connection will become even more intimate after the 
results of the present paper. It was first solved by Kac and van-Moerbeke in [24], 
using the discrete version of inverse scattering of Flaschka. Moser in [29^ solved this 
system explicitly in a fashion similar to |28j . i.e. by construction of action-angle 
coordinates, using the Weyl function and continued fractions. We will not exam- 
ine this result in the paper since it is similar (although more complicated) to the 
corresponding one for the Toda lattice which we describe in Section 2. 

The purpose of this paper is to understand completely the following diagram: 



Toda(g,p) ^ Toda(a,6) 



Invol. 



■ 



Invol. 



Q 

Voltcrra(q) > Volterra(a) 



Dedicated to H. Flaschka and J. Moser. 
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The top arrow is the well-known symplectic realization of the Toda lattice where 
F is the Flaschka transformation. We can think of it as the following symplectic 
bi-hamiltonian realization: 

F : {J\,J 2 ,hx,h 2 ) -> (7ri,7r 2 ,-£/i,-ff 2 ) , 
where tt\ and ir 2 are the well-known linear and quadratic Toda brackets and H\ , H 2 
are the trace of the Lax matrix and the Hamiltonian respectively. They form a bi- 
Hamiltonian pair: nidH 2 = n 2 dHi. Similarly J\ is the standard symplectic bracket 
while J 2 is the Das-Okubo bracket in (q,p) coordinates, hi and h 2 correspond to 
the sum of the momenta and the Hamiltonian respectively. Here, we also have 
a bi-hamiltonian formulation J±dh 2 = J 2 dh±. The multi- hamiltonian structure in 
(a, b) coordinates was obtained in [BJ, [S] while in [q,p) coordinates in [T7j using a 
Theorem of Oevel. There is a recursion operator defined by TZ = J 2 Jf . 

The bottom arrow is a new symplectic bi-hamiltonian realization of the KM- 
system. Let us denote it by 

G : (w 2 ,w 3 ,ii,i ) -> (v 2 ,v 3 , h, I ) . 

The brackets v 2 and V3 are the well-known quadratic and cubic brackets for 
the KM-system defined in [7J. The functions It are defined by I\ — |tr£ 2 , and 
Iq = \og\det(L)\, where L is the Lax matrix. The phase space for the KM-system 
V a is of odd dimension (ai, a 2 , . . . , a 2n +i)- Letting N = 2(n + 1) the dimension of 
V a is N — 1, while the dimension of Tt a< b) is 2N — 1 and the dimension of T/ qtP > is 
2N. 

We consider the space V q in the variables (qi, q 2 , . . . , q^-). On this space we 
define a symplectic Poisson bracket w 2 by the formula 

{quqj} = 1 V i<j . 
This is a constant symplectic bracket and the matrix {qi, qj} has determinant one. 
This bracket w 2 corresponds to the quadratic bracket v 2 via the mapping 

G( qi ,...,q N )=(e^,...,e^-^) . 

We then define a bracket in H N which is mapped to U3 under the transfor- 
mation G. The bracket W3 is defined by the formula: 

{q u qj } = e «-i-« + (1 _ 5 i+1J )e«'-«+ I + e*" 1 "^ + e *-^+ 1 . 
Whenever a term is not defined we omit that whole term. Define 

io = qi - q2 + 93 - 94 H h gjv-i - 9w 

and the Hamiltonian i\ given by 

JV-l 

il=J2 e Ql ~ qi+1 ■ 

i=l 

Then we have the following bi-hamiltonian pair: 

w 2 di\ = w^dio . 

Finally let us comment the two vertical arrows of the diagram. The right hand 
arrow is described in |fOj where the multi-hamiltonian structure of the KM-system 
is obtained by reducing the Toda hierarchy using a Theorem of Fernandes and Van- 
haecke [T8]. The left hand arrow is similarly a reduction of the multi-hamiltonian 
structure of the Toda lattice in (q,p) coordinates to a new multi-hamiltonian struc- 
ture for the KM-system in V q space. We use the involution i/j(q,p) — (q, —p) which 
is a Poisson automorphism of both J 2 and J4 = TZ 2 J 2 . The reduction of J 2 and J4 
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produces the brackets W2 and W3 respectively. We should mention that Moser had 
an algorithm for going in the opposite direction of the right hand side arrow, i.e. 
a procedure for obtaining Toda equations from the Volterra equations by squaring 
the Lax matrix and some chopping. This is described briefly in Section 10. 

2. Moser's solution of the Toda lattice 
The Toda lattice with Hamiltonian function 



N N-l 



N , N-l 

(1) H(q 1 ,...,q N ,p 1 ,...,p N ) 

i=l ~ i=l 

can be transformed via a change of variables due to Flaschka 
of the form L = [B, L], where L is the Jacobi matrix 



to a Lax pair 







fbx 


a 1 


• 


\ 






ai 


b 2 


a 2 ■ 




(2) 


L = 





«2 


h ' 














a-N-i 






u 









(3) 



and B is the projection of L into the skew-symmetric part of L in the Lie algebra 
decomposition, lower triangular plus skew-symmetric. This is an example of an 
isospectral deformation; the entries of L vary over time but the eigenvalues remain 
constant. It follows that the symmetric polynomials of the eigenvalues, Hi = \tr L 1 
are constants of motion. The differential equations in the case of symmetric tri- 
diagonal Lax matrix are: 

cm = ^ (6 l+ i - bi) 
k = 2 (a? - ati) ■ 

Moser's elegant solution of the open Toda lattice uses the Weyl function /(A) and 
an old (19th century) method of Stieltjes which connects the continued fraction of 
/(A) with its partial fraction expansion. The key ingredient is the map which takes 
the (a, b) phase space of tridiagonal Jacobi matrices to a new space of variables 
(Ai, ri) where Ai is an eigenvalue of the Jacobi matrix and r.; is the residue of /(A). 
We present a brief outline of Moser's construction. 

Moser in [28] introduced the resolvent 



R{\) = (XI - L)- 1 , 
and defined the Weyl function 

/(A) - Rnn (A) = (R(X)e N , e N ) , 
where ejv = (Q, 0, . . . , 0, 1). 

The function /(A) has a simple pole at A = Ai and so it admits a partial fraction 
expansion 



N 2 
i—l 

with positive residue r|. It is clear from the formula for calculating the inverse 
of a matrix that /(A) = . We denote by the k by k sub-determinant 
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obtained by deleting the last N — k rows and columns of XI — L. It follows that 
liniA^oo A/(A) = 1 and therefore the residue of /(A) at infinity is — 1. As a result, 

JV 

E^ = i- 

i=l 

More generally, one has a recursion formula: 

(5) A k = (A - &fc)A fc _i - al^Ak-2 ■ 

Moser notes that the mapping (f> between 



T {a,b) = {(oij • • • , a N -i, h, . . . , b N , di > 0)} 

and 

N 

T (\,r) = {(Ai, ...,\N,n,...,rN, Ai < A 2 < ■ • • < Ajv, ^r 2 2 = l, r< > 0)} 

i=l 

is one to one and onto. The inverse mapping cj)^ 1 : Tn r ) ~ * ^~(a,6) corresponds to 
the inverse scattering transform. 

Moser derives the differential equations in the variables (A, r) using the Lax pair 
of the system and the recursion formula ([5|) . The equations take the form 



A, - 

N 

(6) h = -(Ai-^Vfjrj. 

3=1 

If one considers as homogeneous variables then the differential equations ^ 
become linear 

Xi = fi = —XiTi . 

The variables af , 6^ may be expressed as rational functions of Ai and Ti using a 
continued fraction expansion of /(A) which dates back to Stieltjes. Since the com- 
putation of the continued fraction from the partial fraction expansion is a rational 
process the solution is expressed as a rational function of the variables (Ai, ri). 
The procedure is as follows: 

The Rnn element of the resolvent, as defined previously, takes the following 
continued fraction representation: 



(7) /(A) = - 2 

X-b N 2* 

A— 6jf-i — - 



Stieltjes described a procedure that allows one to express a; and bi in terms of 
Ai, . . . , Xn and rj., . . . , r^r. We briefly describe the method. We expand the partial 
fraction expansion of /(A) given in ((4|) in a series of powers of j. 

The coefficient of A 3+1 is denoted by Cj and equals, 

N 

' »< A <- 3 -0,1,.... 

i=l 

The formulas of Stieltjes involve certain i X i determinants which we now define: 
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A, 



CO 
Cl 



Cl 
C2 



Cj-1 
Ci 



Cl 

C2 



C2 
C3 



Ci 
Cj+1 



Cj-1 Ci ... C2i-2 Ci Ci+i . . . C2t_l 

The formulas that give the relation between the variables (a, b) and (r, A) are, 



Aj_iA t+1 

^4i-Bi-2 



Ai-iBi 

N+1 % AiBi_\ 
where A = 1, B = 1, S_i = 0. 
For example, in the case N = 2 



1, 
1,. 



Ai = co, A 2 



CQC2 



and therefore 



Thus, 



a\ = A 2 , bi 



■ 



-4, 



Bi = m, B 2 



ClC 3 



Bo 



B\ A 2 B\ 



b 2 = B l . 



(8) 



r 1 r|(A 2 -A 1 ) 2 
(r? +rf) 2 
A 2 + r|Ai 



r? Ai + r|A 2 



' 1 I" ' 2 

We have written the solution in this form to show that are homogeneous coordi- 
nates, i.e. the solution does not change if we replace r"j by a non-zero multiple. As 
we mentioned earlier in these homogeneous coordinates the differential equations 
become linear. One can check that the differential equations ^ correspond via 
transformation © to the A 2 Toda equations 



di = ai{b 2 - bi) 
t>i = 2a\ 
b 2 = -2a? . 

As Moser notes, it is not too hard to obtain explicit expressions for N = 3 but 
the general case is quite complicated. 

3. Background 

Consider a differential equation on a manifold M defined by a vector field %■ A 
vector field Z is a symmetry of the equation if 

[z, x ] = o. 

A vector field Z is called a master symmetry if 

[[z,x],x] = o, 

but 

Master symmetries were first introduced by Fokas and Fuchssteiner in |2f j in con- 
nection with the Benjamin-Ono Equation. 
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A bi-Hamiltonian system is defined by specifying two Hamiltonian functions H\ , 
H 2 and two Poisson tensors tt± and -k 2 , that give rise to the same Hamiltonian 
equations. Namely, ttiV H 2 = ^Vffi. The notion of bi-Hamiltonian structures is 
due to Magri [26]. Suppose that we have a bi-Hamiltonian system defined by the 
Poisson tensors 7Ti, ~k 2 and the Hamiltonians Hi, H 2 . Assume that 7Ti is symplectic. 
We define the recursion operator 1Z = 7r 2 TTi , the higher flows 

Xl = ft*" 1 *! , 

and the higher order Poisson tensors 

7T- = K^tti . 

For a non-degenerate bi-Hamiltonian system, master symmetries can be generated 
using a method due to Oevel [3"U] , 

Theorem 1. Suppose that Xq is a conformal symmetry for both tti, ir 2 and Hi, 
i.e. for some scalars \, fj,, and v we have 

C-XqKi = A-7Ti, Cx ^2 — ^2, Cx Hi = vH\. 

Then the vector fields AQ — TVXq are master symmetries and we have, 

(a) C Xi Hj = (u+U-l + i)^- X))H i+j 

(b) C Xi nj = (/* + (j - i - 2)(/x - A))7r l+i 

(c) [XuXj] - (»-X)(j-i)X i+j . 

A symplectic realization (see [37]) of a Poisson manifold (M, tt) is a symplectic 
manifold (S, oj) together with a surjective Poisson submersion / : S — > M. In this 
paper our realizations will have additional structure; they will be symplectic bi- 
hamiltonian realizations as in [33] . [34] . Suppose on M we have a bi-hamiltonian 
pair 7Ti, 7T2 and two Hamiltonians H±, Hi such that nidHi = itidHi. A symplectic 
bi-hamiltonian realization will be a manifold S, of even dimension, together with 
two Poisson tensors J±, J2 with J± symplectic and a surjective submersion F : S — > 
M which is a Poisson mapping between and 7T,, i.e. {F* f, F*g}j t — F*{f,g}^ i 
for all f,g£ C°°(M). In addition, Jidhi — Jidhi where hi = HiO F. 

We shall see that the relation between the Toda and Volterra systems relies 
on special symmetries of the phase spaces. We will need a theorem of Fernandcs 
and Vanhaecke which gives conditions under which the fixed point set of a Poisson 
action inherits a Poisson bracket. 

Although we will be interested mainly in finite symmetries, we have the following 
general result [TB] : 

Theorem 2. Suppose that (Af, {•,■}) is a Poisson manifold, and G is a compact 
group acting on M by Poisson automorphisms. Let N — M G be the submanifold of 
M consisting of the fixed points of the action and let it, : N <—* M be the inclusion. 
Then N carries a (unique) Poisson structure {•, -} N such that 

(9) i*{F u F 2 } = {i*Fi,i*F 2 } N 

for all G-invariant functions F±, F 2 € C°°(M). 

Remark 3. The previous result can be seen as a particular case of Dirac reduction 
(for the general theorem on Dirac reduction, see Weinstein f[37j . Prop. 1.4) and 
Courant ([5\, Thm. 3.2.1). 

Since this result applies in particular when G is a finite group, we have: 

Corollary 4. Suppose that (Af, {•,•}) is a Poisson manifold, and G is a finite group 
acting on M by Poisson automorphisms. Then the fixed point set N — M G carries 
a (unique) Poisson structure {•, -} N satisfying HJ). 
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Let us consider the special case G — Z2. Then G — {I, <fi}, where tj> : M — * M 
is a Poisson involution. We conclude that N = M G = {x : 4>(x) = x} has a unique 
Poisson bracket satisfying equation ([9]). So we see that Theorem [2] contains as a 
special case the following result, which is known as the Poisson involution theorem 
(see [HI [H]). 

Corollary 5. Suppose that (M, {•,•}) is a Poisson manifold, and <fi : M — > M is 
a Poisson involution. Then the fixed point set N — {x E M : (f>(x) — x} carries a 
(unique) Poisson structure {•, -} N such that 

i*{F 1 ,F 2 } = {i*F 1 ,t*F 2 } N 

for all functions F±,F2 € C°°(M) invariant under <f>. 



4. KM-SYSTEM 

The Volterra system, also known as KM system is defined by 



(10) 







i = 1,2, 



where clq — a n +i =0. It was studied originally by Volterra in [35j to describe 
population evolution in a hierarchical system of competing species. In [29] Moser 
gave a solution of the system using the method of continued fractions and in the 
process he constructed action-angle coordinates. Equations (fT0|) can be considered 
as a finite-dimensional approximation of the Korteweg-de Vries (KdV) equation. 
They also appear in the discretization of conformal field theory; the Poisson bracket 
for this system can be thought as a lattice generalization of the Virasoro algebra 
|16j . The variables at are an intermediate step in the construction of the action- 
angle variables for the Liouville model on the lattice. 

The Volterra system is usually associated with a simple Lie algebra of type A n . 
Bogoyavlensky generalized this system for each simple Lie algebra and showed that 
the corresponding systems are also integrable. See [H [3] for more details. The 
relation between Volterra and Toda systems is also examined in [TU], [TT] . 

The Hamiltonian description of system (|T0|) can be found in ,15] and . The 
Lax pair is given by 







L 


= [B,L], 




where 












(° 


1 


••• 


••• ON 




ai 





1 




(11) L = 





a 2 


'•• 







Vo 




... 


'■■ 1 

a n 0/ 



and 
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/ 





B = 



-a\a,2 
: -a 2 a 3 



^TrL 2 * are constants of motion. The system 



V -a n _ia„ 0/ 

It follows that the functions = ^Tri 21 
pop is integrable for any value of n but in this paper we restrict our attention only 
to the case where n is odd. 

Following [7J we define the following quadratic Poisson bracket, 

{di, a,i + i} — didi+i, 

and all other brackets equal to zero. We denote this bracket by V2- For this bracket 
detL is a Casimir and the eigenvalues of L are in involution. Of course, the functions 
Ii are also in involution. Taking the function l\ = di as the Hamiltonian we 
obtain equations (|T0|) . i.e. V2d.l1 is the KM Hamiltonian vector field. 

In [7] one also finds a cubic Poisson bracket which corresponds to the second 
KdV bracket in the continuum limit. It is defined by the formulas, 

Wi, d i+ i} = dia i+ i{d t + a i+ i) 

{di,di+2} = 0,idi + idi + 2 , 

all other brackets are zero. We denote this bracket by V3. In this bracket we still 
have involution of invariants. We also have Lenard type relations of the form 

«3 dhi = V 2 dI 2i+ 2 ■ 

In [7] appears a bracket that is homogeneous of degree one, a rational bracket 
constructed using a master symmetry. This bracket, denoted by vi, has 1\ = ^TrL 2 
as Casimir and the Hamiltonian is I2 = jTiL 4 . The definition of the bracket is the 
following: We define the master symmetry Y_i to be 



where the fi are determined recursively as follows, 



fi 



-1, h 



«2 



02i-l 



—f2i-l, fli-\ 



1. 



Taking the Lie derivative of V2 in the direction of K_i we obtain v±, a Poisson 
bracket that is homogeneous of degree 1. For n = 5, v\ takes the form: 



(12) 



{ax, a 2 } 


= a-2 


{ai,a 3 } = 


-02 


{^2,03} 


= a 2 


{o 2 , a 4 } = 


a2^4 


{a 3 ,a 4 } 


= d4 


{a 3 ,a 5 } = 


— 04 



{ax,a 4 } = ^ 



{a2,«5} 
{a 4 , a 5 } : 



Note that the KM Hamiltonian can be expressed as V2dli = and this is 

a bi-hamiltonian formulation of the KM flow. In this paper we rediscover these 
brackets using a recursion operator. 

The higher Poisson brackets are constructed in [7J using a sequence of master 
symmetries Yi. For example, the bracket V2 is obtained from v\ by taking the 
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Lie derivative in the direction of the first master symmetry Y\ . Similarly, the Lie 
derivative of V2 in the direction of Y\ gives V3 . 

The brackets v\, V2 and V3 are just the beginning of an infinite hierarchy which 
we will derive again in this paper using a different method. 

5. Toda Lattice 

In this paper we will deal with the finite, non-periodic version of the Toda lattice. 
However, there are also other interesting variations worth mentioning. There is a 
generalization due to Deift, Li, Nanda and Tomei fJJ] who showed that the system 
remains integrable when L is replaced by a full (generic) symmetric n x n matrix. 
The functions Hi = i Tr L % are still in involution but they are not enough to 
ensure integrability. In other words, the existence of a Lax pair does not guarantee 
integrability. There are, however, additional integrals which are rational functions 
of the entries of L. The method used to obtain these additional integrals is called 
chopping and was used originally in [14] 

The classical Toda lattice was also generalized in another direction. One can 
define a Toda type system for each simple Lie algebra. The finite, non-periodic 
Toda lattice corresponds to a root system of type A n . This generalization is due 
to Bogoyavlensky [3]. These systems were studied extensively in [25] where the 
solution of the system was connected intimately with the representation theory of 
simple Lie groups. There are also studies by Olshanetsky and Perelomov [31] and 
Adler, van Moerbeke [T]. The description of the systems, following p] and [25J is 
as follows: 

Let g be any semi-simple Lie algebra, equipped with its Killing form (• | •). One 
chooses a Cartan subalgebra t) of g, a root system A = A(t), g) of f) in Q, a basis of 
simple roots LI, and a set of positive roots A+. 

The Lax pair (L(t), B(t)) in g can be described in terms of the root system as 
follows: 

1 1 
L(t) = J2 b i(t)H ai +J2ai(t)( x *i + x -*i) , 

I 

i=l 

where H ai is an element of f) and X ai is a root vector corresponding to the 
simple root a,. 

We consider now the classical Toda lattice but we will use a Lax pair slightly 
different than f2}. 

Let D be the diagonal matrix with entries di where d\ = 1 and di = a\a2 ■ ■ ■ a^_i i — 
2, 3, ... TV. In [55J Kostant conjugates the matrix L in @ by the matrix D and the 
resulting matrix DLD^ 1 has the form 



(bi 1 





ai 6 2 


1 


a 2 


63 



: '•• '■• '•■ 

: '■• '■• 1 

\0 ajv-i b N J 

Denoting this matrix again by L the equations take the Lax form 



10 



PANTELIS. A. DAMIANOU 



L(t) = [L(t),PL(t)] 

where P is the projection onto the strictly lower triangular part of L(t). The de- 
composition here is strictly lower plus upper triangular. This form is convenient in 
applying Lie theoretic techniques to describe the system. Note that the diagonal 
elements correspond to the Cartan subalgebra while the subdiagonal elements cor- 
respond to the set II of simple roots. The equations of motion for the Toda lattice 
(in Kostant form) are: 



a>i = a,i (bi+i - hi) 

The functions H % = iTr U are independent invariants in involution. 
There exists a Lie-Poisson bracket given by the formula 

(15) , = ~ ai 

{ai,b l+ i\ = <n ; 

all other brackets are zero. H\ = b\ + 62 + ■ ■ • + b^ is the only Casimir. The 
Hamiltonian in this bracket is H2 — \ tr L 2 . We also have involution of invariants, 
{Hi, Hj} = 0. The Lie algebraic interpretation of this bracket can be found in [25] , 
We denote this bracket by -k\. 

The quadratic Toda bracket appears in conjunction with isospectral deformations 
of Jacobi matrices. It is a Poisson bracket in which the Hamiltonian vector field 
generated by H\ is the same as the Hamiltonian vector field generated by H% with 
respect to the tti bracket. The defining relations are 



(16) 



{a i7 a i+ i} = a t a i+1 

{di,bi} = —dibi 

{di,b i+1 } = dib i+ i 

{biM+i} = a i ; 



all other brackets are zero. The bracket TT2 is easily defined by taking the Lie deriv- 
ative of tv 1 in the direction of suitable master symmetry X\, see [S] for details. This 
bracket has det L as Casimir and H\ = tr L is the Hamiltonian. The eigenvalues of 
L are still in involution. Furthermore, tt2 is compatible with tt\. We also have 

■KidHi = TTidH l+1 . 

Finally, we remark that taking the derivative of tt2 in the direction of X\ yields 
another Poisson bracket, 7r3, which is cubic in the coordinates. The defining rela- 
tions for 7r 3 are 



(17) 



{di,di+i} = 2diCLi+ibi+i 

{di,bi} = -dib\-d\ 

{a,i,b i+ i = aibf +1 +a 2 

{di,bi + 2} = didi+i 

{a>i+i> bi} = —dicii+i 

{bi,b i+ i} = di (bi + b i+ i) ; 



all other brackets are zero. The bracket 113 is compatible with both 7Ti and 7T2 and 
the eigenvalues of L are still in involution. The Casimir for this bracket is tr_L _1 . 

In fact there is an infinite hierarchy of Poisson tensors 71^ , master symmetries Xi 
and invariants Hi and they obey some deformation relations. We quote the results 
from refs. 0, 0. 
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Theorem 6. 

i) Kj, j > 1 are all Poisson. 

ii) The functions Hi, i > 1 are in involution with respect to all of the TTj. 
in) Xi (Hj ) = (i + , t > -1, j > 1. 

roj X^ttj = (i - i - 2)7r j+J -, i > -1, i > 1. 

«; [x^ x,-] = (j - i)x i+jl i > o, j > o. 

vi) TTj dHi = TTj—i dHi + \. 

6. TODA LATTICE IN (q,p) COORDINATES 

Let Ji be the symplectic bracket with Poisson matrix 



P 
-I 0y 

where / is the N x N identity matrix. The bracket J\ is mapped precisely onto 
the bracket tti under the Flaschka transformation F: 

(18) a i = e*-*+ 1 , bi = - Pi . 

We define the following tensor, J2, due to Das and Okubo 13J : 



(19) J 2 = 



where A is the skew-symmetric matrix defined by = 1 = — &ji for i < j, B is the 
diagonal matrix (— pi, — P2, ■ ■ ■ , ~Pn) and C is the skew-symmetric matrix whose 
non-zero terms are c^+i = — Cj+i^ = e 9i ~ 9i+1 for i = 1, 2, . . . , N — 1. The bracket 
J2 is mapped precisely onto the bracket 712 under the Flaschka transformation (|18|) . 
It is easy to see that we have a bi-Hamiltonian pair. We define 

hi = -(pi +P2 H h Pa) , 

and /i2 to be the Hamiltonian: 

N jV-1 



i=l i=l 

Then we obtain the bi-hamiltonian pair 



Jidh2 = Jidhi . 
We define the recursion operator as follows: 

n = J2J1 1 • 

The matrix form of 1Z is quite simple: 



(20) 



Using the recursion operator we obtain the higher order Poisson tensors 

Ji =n i ~ 1 Ji i = 2,3, . . . . 
Following }17j we define the conformal symmetry 

Ti: 

dpi 



N 1 N 1 

Z = J2(N-2i + l)— + ^ 1 



i—l 

It is straightforward to verify that 

C-z Ji = — Ji 
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£z J 2 = . 

In addition, 

Z (hi) = hi 
Z (h 2 ) = 2h 2 . 

Consequently, Zq is a conformal symmetry for Ji, J 2 and h\. The constants 
appearing in Theorem Q] are A = — 1, /i = and i/ = 1. According to Oevel's 
Theorem we end up with the following deformation relations: 

[Z h hj] = (i + j)h i+j 
LziJj = (i — i — 2)Ji+j 
= (j -i)Z i+3 ■ 

Switching to Flaschka coordinates, we obtain relations hi)- v) of Theorem H>] 

7. From Toda to Volterra 

We consider T( a ,b) the phase space of the Toda lattice in Flaschka coordinates 
and the space V a of KM-system in a coordinates. Note that V a is not a Poisson 
subspace of T^ a b y However, as was demonstrated in |10) V a is the fixed manifold 
of the involution <f> : Tr a ,b) ~ > T( a ,b) defined by 

<f)(ai, a 2 • • ■ , ajv-i,&i,&2 • • • , &jv) ^ (ai,a 2 . . ., a,N-i,—bi, -b 2 . . . , -b N ), 

and we have the following result: 

Theorem 7. : Tr a ,b) ~ * ^(a,fc) is a Poisson automorphism of (Tr a irk), if k is 
even. 

Therefore, by Corollary El V a inherits a family of Poisson brackets 1)2,1)3,.... For 
example, the quadratic bracket can be computed from formulas (TTo| and is given 
by 

(21) {ai,aj} = ti.o , . 1 - 5 i+ ij), 

while the Poisson bracket V3 is found to be given by the formulas 

(22) {ai,cn + i} = aia i+ i(ai + a,i + i), (i = 1, . . . , n - 1) 
{ai, a l+2 \ = aia i+ ia l+2l (i = 1, . . . , n - 2) 

i.e. we obtain v 2 and U3 of Section 4. It follows that the restriction of the 
integrals H 2 k to V a gives a set of commuting integrals, with respect to these Poisson 
brackets. Also, the Lax equations (fT5)) lead to Lax equations for the corresponding 
flows, merely by putting all bi equal to zero. We recover the vector field 

(23) d,- = Oj(aj_i — eti+i), i = l,...,n, 

and a family of integrable systems admitting a multiple hamiltonian formulation: 

v k dli = Vk-idli+i, (k = 1,2, ...), 
where Vk is the restriction of n 2 k to V a and // is the restriction of H 2 \ to V a . 
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8. Symplectic realization 

Let N = 2(n + 2) and consider coordinates (qi, (72, • ■ ■ , <?7v) in R N ■ We define the 
following transformation G from R 2 ™+ 2 to R 2n+1 , 

(24) a i = e qi - q *+ 1 i = 1, 2, . . . ,N - 1 . 

The Hamiltonian in q coordinates is given by 



N-l 

(25) h = e qi ~ qi + x . 

i=i 

We define the Poisson bracket 

Ui,Qj} = 1 V i<j . 

Let us denote this constant Poisson tensor by W2 (since it is the analogue of V2 in 
q coordinates). 

Hamilton's equation W2<iii become 

(26) = _ e «-i-« _ e »-»+i . 

It is straightforward to check that Hamilton's equations correspond in the a— space 
to the KM-systcm {TDJ via the mapping G : R N -> R^ 1 



(27) G( qi) ...,q N )=(e q i- q *,...,e q »-i- q ») . 

In fact we calculate: 

d i = e*-*+ 1 (<h-g i+ i) 

= aj(aj+i — a^—x) ■ 

The symplectic bracket u>2 in space corresponds to the quadratic bracket V2 
in V a space. 

We then define a bracket W3 in which is mapped to W3 under the transfor- 
mation G. The bracket W3 is defined by the formula: 

qj } = e q *-i- q * + (1 - 5 l+lj )e q '- q '+ 1 + tfh-i-ii + e «i-*+ 1 . 

Whenever a term is not defined we omit that whole term. This happens only when 
i = 1 or j = N. We note that W2 is compatible with W3. 
Define 

N 

k = *^2{-i) k+1 qk = qi - <72 + q-3 - qi H 1- <in-i - qn ■ 

k=l 

Remark 8. It is not difficult to discover iq. According to |12j . in the presence of an 
invertible Nijenhuis tensor a natural choice of functions to form a bi-hamiltonian 
pair is 

^log(det7?.) and ^txTZ. It turns out that the determinant of the Toda recursion 
operator restricted to the Volterra phase space is e 2l ° and the trace equals 2i\. 

Then 

^{g.,ioU = -e'- 1 -''-e« +1 . 
In other words we have a bi-hamiltonian pair 
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Wc define a recursion operator as follows: 

TZ = W3W2 1 ■ 

In q coordinates, the symbol Xj is a shorthand for xu ■ It is generated, as usual, by 

Xi = K i - 1 X i. 

Note that i\ corresponds under mapping (|27|) to a constant multiple of I\ = 
^Tr (L) 2 . In a similar fashion we obtain the higher order Poisson tensors 

w,i = TV i = 3, 4, 

We finally define the conformal symmetry 

N „ 

X =J2(N-i + l)- 

i=l Hl 

The Poisson tensors W2, W3 and the functions io, %\ define a bi-Hamiltonian pair. It 
is straightforward to verify that 

£x W2 = 0, C x „w 3 = w 3 , £ Xo ii = h- 
Consequently, Xq is a conformal symmetry for W2, w 3 and i%. The constants ap- 
pearing in Oevel's Theorem are A = 0, jjL = 1 and v = 1. Therefore, we end up with 
the following deformation relations: 

[Xk,ij] = (k + j)ik+j 

L x k Vj = (j - k — 2)v k+j 

[X k ,X J ] = (j-k)X k+] . 
Projecting to the a— space under mapping (|27p we obtain the multiple hamiltonian 
structures of [7j. 

Remark 9. Note that we may define a Poisson tensor wi by the formula wi — 
W2W 3 ~ 1 ui2 as in]§\. The projection ofw\ under transformation \2T\j gives v% of Sec- 
tion 4- Perhaps a more appropriate bottom arrow in the diagram of the introduction 
is 

G : (l01,1i>2,*l,i2) -» (vi,V2,Ii,I 2 ) ■ 

9. From Toda (q,p) space to Volterra q space 

In this Section we will explain the origin of the symplectic bi-hamiltonian real- 
ization of Section 8. The idea is to use a Poisson involution in the Toda (q,p) space 
and to reduce the equations to the Volterra V q space. We consider TV 9)P ), the phase 
space of the Toda lattice in natural coordinates. In Section 6 we have constructed 
a bi-hamiltonian system given by the Poisson tensors J\, J2 and the Hamiltonians 
h\ and /i2- We also have a recursion operator TZ which gives rise to a sequence of 
Poisson tensors Ji i = 1, 2, 

Define the involution ip : Tr q , p ) — * T(q,p) by the formula 

if}(q 1 ,...,q N ,p 1 ,...,p N ) = {qi, . ■ ■ ,qN,-Pl, • • • ,-Pn) ■ 
We have the following result: 

Theorem 10. ip : Tr qp -\ — > Tr qp -\ is a Poisson automorphism of (T/ q p j , J2k), k — 
1,2,.... 
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We will not give the proof of this result since it is entirely analogous to the proof 
of Theorem [7] which is given in |10j . 

Therefore using Corollary [5j V q inherits a family of Poisson tensors W2, ti>3, 

For example, the bracket u>2 is clearly the A block of the Poisson matrix J2 in (fT9|) . 
On the other hand it is straightforward to compute that in the bracket J4 we have 



Qj] = pf + P l Pj + Pj + e^ 1 9i + (1 - S i+ i,j)e' 
Therefore the reduction of J4 to V q is precisely W3 of Section 8. 



+ e 



10. Moser's recipe 

Moser in [29] describes a relation between the KM system (flQ|) and the non- 
periodic Toda lattice. The procedure is the following: Form L 2 which is not any- 
more a tridiagonal matrix but is similar to one. Let {ei, e2, . . . , e n } be the standard 
basis of R™, and E a = {spane2i_i, i = 1,2,...}, E e = {spane2i, i = 1,2,...}. 
Then L 2 leaves E a , E e invariant and reduces to each of these spaces to a tridiag- 
onal symmetric Jacobi matrix. For example, if we omit all even columns and all 
even rows we obtain a tridiagonal Jacobi matrix and the entries of this new matrix 
define the transformation from the KM-system to the Toda lattice. We illustrate 
with a simple example where n — 5. 

We use the symmetric version of the KM system Lax pair (due to Moser) given 

by 



(28) 



L = 



/0 

a 1 





01 


«2 








«2 











«4 








It is simple to calculate that L 2 is the matrix 



/ a 2 





aia 2 





\ 





a\ + al 





a 2 a 3 





aia 2 





a\ + a\ 





a 3 a 4 





0203 





a 2 ^ + a\ 





V 





a 3 a 4 





< ) 



(29) 



Omitting even columns and even rows of L 2 we obtain the matrix 

(a\ a\a2 
0102 a\ + al a 3 az 
a 3 a 4 a\ 

This is a tridiagonal Jacobi matrix. It is natural to define new variables A\ — a\a2, 
A2 = 0304, B\ — a\ , B2 = al + a 2 , B3 = a 2 . The new variables A\, A2, B\, B2, B3 
satisfy the Toda lattice equations ([3|). 

This procedure shows that the KM-system and the Toda lattice are closely re- 
lated: The explicit transformation which is due to Henon maps one system to the 
other. The mapping in the general case is given by 



(31) 



1 



Aj — — — y/a,2ia,2i-i 



0-2i 



The equations satisfied by the new variables A4, Bi are given by: 



A t = Ai (B i+1 - Bi) 
B t = 2{A 2 -AU). 
These are precisely the Toda equations in Flaschka's form [19 . 
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This idea of Moser was applied with success to establish transformations from the 
generalized Volterra lattices of Bogoyavlensky [2 [3] to generalized Toda systems. 
The relation between the Volterra systems of type B n and C n and the corresponding 
Toda systems is in [10] . The similar construction of the Volterra lattice of type 
D n and the generalized Toda lattice of type D n is in [TT] . It turns out that the 
Volterra D n system corresponds not to the Toda D n system but to a special case 
of the Sklyanin lattice [32] . 
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